THE EYNARD-ORANTIN RECURSION FOR THE TOTAL ANCESTOR 

POTENTIAL 



TODOR MILANOV 

Abstract. It was proved recently that the correlation functions of a semi-simple 
cohomological field theory satisfy the so called Eynard-Orantin topological re- 
cursion. We prove that in the settings of singularity theory, the relations can be 
expressed in terms of periods integrals and the so called phase forms. In particular, 
we prove that the Eynard-Orantin recursion is equivalent to N copies of Virasoro 
constraints for the ancestor potential, which follow easily from the definition of 
the potential. 
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1. Introduction 

The Eynard-Orantin recursion (see (51) was discovered first for the correlation 
functions of certain matrix integrals. However, its applications go beyond the theory 
of matrix models. The recursion is turning into a powerful tool for computing the 
correlation functions in various quantum field theories. In particular, it provides 
an efficient algorithm for computing quite complicated invariants such as Gromov- 
Witten invariants and certain polynomial invariants in knot theory. 

In order to set up the recursion one needs an analytic curve, called spectral curve, 
two holomorphic functions on it, and a certain symmetric 2-form satisfying some 
additional properties. At first, one might think that this is a serious restriction, so 
the applications would be only limited. The surprising fact however, is that the ini- 
tial data can be determined from the 1 -point and the 2-point correlation functions 
only (of a given quantum field theory) (see 0). This feature makes the recursion 
quite universal. In particular, this observation was exploited in the paper fl4], where 
the authors prove that the ancestor Gromov-Witten (GW for shortness) invariants 

l 
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of manifolds with a semi- simple quantum cohomology can be computed via the 
Eynard-Orantin recursion. Although our work appears after Q, the main observa- 
tion namely, that one should study the n-point series © and that they should satisfy 
the Eynard-Orantin recursion with kernel given by formulas (|36l ) and (|37l) was done 
independently. 

The goal in this paper is to interpret the Eynard-Orantin recursion in terms of 
differential operator constraints for the total ancestor potential. In particular, this 
allows us to obtain a simple proof of the recursion relation. In particular, we prove 
that the correlation functions can be expressed in terms of period integrals and phase 
forms, which suggests that they should be compared to the correlation functions of 
the twisted Vertex algebra representation introduced in [0. 

1.1. Preliminary notation. Let / e O c im be the germ of a holomorphic function 
with an isolated critical point at 0. We fix a miniversal deformation F(t, x),t e B and 
a primitive form to in the sense of K. Saito lfT3l[T5l . so that B inherits a Frobenius 
structure (see IfTTl O). In particular we have the following identifications (c.f. 
Section O): 

T*B = TB = Bx T Q B = B x H, 

where H is the Jacobi algebra of /, the first isomorphism is given by the residue 
pairing, the second by the flat residue metric and the last one is the Kodaira-Spencer 
isomorphism 

T B = H, d/dti » d ti F[ =0 mod (/„,..., f m ). (1) 
We need also the period integrals 

lf(t,A) = -d B (InY'dY f d~ l a) 6 T;B = H, (2) 

where a is a cycle from the vanishing cohomology, d B is the de Rham differential 
on B, and d~ l a> is any (n - l)-form 77 such that drj = co. The periods are multivalued 
analytic functions on B x P 1 with poles along the so called discriminant locus (c.f. 
Section [2^21) . We make use of the following series 

P(f, A; z) = J] ^(t, A) {-z)\ </> a {t, A;z) = J] V, A)(-z)* dA. 

keZ keZ 

Note that <ff(t, A; z) = d pl f a (t, A; z). 

Let B ss c B be the subset of semi-simple points, i.e., points t e B such that the 
critical values of F(t, •) form a coordinate system in a neighborhood of t. For every 
t e B ss Givental's higher-genus reconstruction formalism gives rise to ancestor 
correlation functions of the following form 

{v^ v n xl/") g , n {t), Vi 6 H, k t e Z + (l < i < n). (3) 
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Apriory, each correlator depends analytically on( e B ss , but it might have poles 
along the divisor B \ B ss . Given n vanishing cycles a\, . . . , a n and a generic point 
t 6 B we define the following n-point symmetric forms 

(D a g ) r a »{t-,A x ,...,A n ) = (f + \t,A x ,ilfd,...A7(t,A n ,iJf n )) g {t), (4) 

where the + means truncation of the terms in the series with negative powers of z. 
The functions © will be called n-point series of genus g. They should be inter- 
preted formally via their Laurent series expansions at the singular points. They are 
the main object of our study and we expect that they have many remarkable proper- 
ties yet to be discovered. Probably the first question to be addressed is whether the 
n-point functions are global objects, i.e., multivalued analytic functions with finite 
order poles on the configuration space CnCP 1 ) = ^(P 1 )/^, where 

^(P 1 ) = {Oli, . . . , A n ) e (P 1 )" : At ± Aj for i ± j] 

and S„ is the symmetric group acting by permutation of the coordinates. 

1.2. Statement of the results. Let us assume that t e B ss is generic so that the 
function F(t, •) has N := dim c H pairwise different critical values Uj(f). Let Bj be 
a cycle vanishing over A = uj. We introduce the following quadratic differential 
operator 

Y% =: (d A &(t, A)) 2 : + A), (5) 

where : : is the normal ordering, the differential operator (c.f. Section 13711) 

d A &(t, A) = d A fj(f, A) + d A f\t, A) 

is defined by 

CO N r, 

d A m, x) + = Y,Yu ( ~ 1)k+l C 1 ^ ^ nm ti • < 6) 

co N 

d x tf>J(!, A)_ = J] J] ffi*' V> v ') r V2 4k . (7) 

k=0 i=l 

where {v'} and {v,} are dual bases for H with respect to the residue pairing (, ). 
Finally, F\ 6 j J '^ J is the free term in the Laurent series expansion of the propagator 

oo 

[d A f + \t,p), dfj(t, A)] = 2(jx- A)' 2 + J] $ 0" " X) k - 

k=0 

The definition © is very natural from the point of view of vertex algebras (see 
lO). It is the field that determines a representation of the Virasoro vertex operator 
algebra with central charge 1 on the Fock space 

C n [[q ,qi + l,q 2 ,...]l q k = {q\, . . . ,q N k ), C,:=C((V^)), 
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where 1 is the unity of the local algebra H. We may assume that Vi = 1. 

Let us denote by 3l t (%; q) the total ancestor potential of the singularity. By defi- 
nition, it is a vector in the Fock space of the form (c.f. Section IX2l) 



where t(i/0 = Tjkj t k v '^ /k anc ^ me re l a ti° n between the set of formal variables {q' k } 
and {f k } is given by the dilaton shift 



Note that although the periods are multi-valued analytic functions, the above ex- 
pression is single valued with respect to the local monodromy around A = Uj, so the 
residue is well defined. Our first result is the following . 

Theorem 1.1. The total ancestor potential satisfies the following constraints: 

L m -ij M t (h; q) = 0, m > 0, 1 < i < N. 

The proof follows easily from the definition of the ancestor potential and some 
known properties of the periods ©• More precisely one can prove that Y" J A &l t is 
regular near A = uj, so each residue vanishes. The main property of the construc- 
tion © is that Lm-\j, has only one term that involves the dilaton- shifted variable q\ 
and this term is q\d/dq' m . This fact allows us to interpret the differential operator 
constraints as a system of recursion relations. Our next result is the following. 

Theorem 1.2. The differential operator constraints determine a system of recursion 
relations that coincide with the local Eynard-Orantin recursion. 

We postpone the definition of the Eynard-Orantin recursion until Section HI Fol- 
lowing @| we give the definition of the recursion only locally. It will be very impor- 
tant to find the global formulation (see 0) as well, but for now this seems to be a 
very challenging problem, except may be for simple and simple elliptic singularities 
in singularity theory or the projective line and its orbifold versions in GW theory. 
The problems is that in the settings of singularity theory, or GW theory, the spectral 
curve has highly transcendental nature. Its description requires inverting the pe- 
riod map: very classical and very difficult problem. Our set up is slightly different 
from the standard conventions ((H 0] [/J 0), because we would like to work with 
multivalued correlation functions. The main idea is that whatever is the spectral 




g,n=0 




(8) 




m > 0, 1 < i < N. (9) 
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curve E, we always have a projection S —* P 1 which in general is an infinite sheet 
branched covering. Galois theory tells us that studying the field of meromorphic 
functions on S is the same as the field of multivalued meromorphic functions on P 1 
invariant under some monodromy (Galois) group. Our proposal is to formulate the 
Eynard-Orantin recursion for correlation functions on P 1 that take values in some 
local system X. The definition that we give in Section @] is just a first attempt to 
set up this idea. Probably a better formulation is possible if one takes into account 
more examples not only the ones that come from GW theory. 

Let us point out that although we work in the settings of singularity theory, the 
differential operators © can be defined for any semi-simple Frobenius manifold 
that has an Euler vector field. The period vectors should be introduced as the solu- 
tions to a system of differential equations (see Lemma |2~3l ) and it is easy to see that 
the proofs of Theorems 11.11 and 11.21 remain the same. In particular, we obtain the 
main result of flU 

Corollary 1.3. The n-point series of a semi-simple cohomological field theory sat- 
isfy the local Eynard-Orantin recursion relations. 

1.3. Acknowledgement. I am thankful to M. Mulase for teaching me how the 
Eynard-Orantin recursion can be constructed in terms of 1- and 2-point functions. 
Also, I thank S. Shadrin for useful e-mail communication and clarifying some 
points from [4J. This work is supported by Grant-In- Aid and by the World Pre- 
mier International Research Center Initiative (WPI Initiative), MEXT, Japan. 

2. Frobenius structures in singularity theory 

Let / : (C 2/+1 , 0) — > (C, 0) be the germ of a holomorphic function with an isolated 
critical point of multiplicity N. Denote by 

H = C[[x ,...,x 21 ]]/(d x J,...,d X2l f) 

the local algebra of the critical point; then dim H = N. 

Definition 2.1. A miniversal deformation off is a germ of a holomorphic function 
F : (C N x C 2/+1 ,0) — > (C, 0) satisfying the following two properties: 

(1) F is a deformation off, i.e., F(0,x) = f(x). 

(2) The partial derivatives dF/dt 1 (1 < i < N) project to a basis in the local 
algebra 

0&,q[[xq, x 2 i]]/(d X0 F, d X2l F). 

Here we denote by t = (t 1 , . . . , t N ) and x = (x , . . . , x 2 i) the standard coordinates 
on C N and C 2?+1 respectively, and O is the algebra of germs at of holomorphic 
functions on C N . 
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We fix a representative of the holomorphic germ F, which we denote again by F, 
with a domain X constructed as follows. Let 

B 2M c C 2/ + 1 ^ g = B N c C N ^ ^1 c c 

be balls with centers at and radii p, 77, and 6, respectively. We set 

S =BxB l s cC N xC, X = (Bx B 2 p M ) n cC^x C 2,+1 , 

where 

cp: Bx B 2 J +l -» fix C, (f,x) i-> (t,F(t,x)). 

This map induces a map : X — > 5 and we denote by X s or the fiber 

^ = *a = {(t, x)eX\ F(t, x) = A}, s = (t,A)£S. 

The number p is chosen so small that for all r, < r < p, the fiber X intersects 
transversely the boundary dB 2l+l of the ball with radius r. Then we choose the num- 
bers n and 8 small enough so that for all s e S the fiber X s intersects transversely the 
boundary dB 2l+1 . Finally, we can assume without loss of generality that the critical 
values of F are contained in a disk 5^ with radius 60 < 1 < 6. 

Let Z be the discriminant of the map 0, i.e., the set of all points s e S such that 
the fiber X s is singular. Put 

S'=S\IcC"xC, X' = f'(S')cXc C*xC 2,+1 . 

Then the map <f>: X' — » 5' is a smooth fibration, called the Milnor fibration. In 
particular, all smooth fibers are diffeomorphic to X l . The middle homology group 
of the smooth fiber, equipped with the bilinear form (-|-) equal to (-1)' times the 
intersection form, is known as the Milnor lattice Q = H 2 i(X y,Z). 

For a generic point s e E, the singularity of the fiber X s is Morse. Thus, every 
choice of a path from (0, 1) to s avoiding E leads to a group homomorphism Q — > 
H 2 i(X s ;Z). The kernel of this homomorphism is a free Z-module of rank 1. A 
generator a e Q of the kernel is called a vanishing cycle if (alar) = 2. 

2.1. Frobenius structure. Let T B be the sheaf of holomorphic vector fields on B. 
Condition (2) in Definition 12 . 1 1 implies that the map 

d/df h-> dF/df mod (d XQ F, d X2l F) (1 < i < N) 

induces an isomorphism between T B and p*O c , where p : X — » B is the natural 
projection (t, x) h-> t and 

O c :=O x /(d X0 F,...,d X2l F) 

is the structure sheaf of the critical set of F. In particular, since Oc is an algebra, the 
sheaf 7~b is equipped with an associative commutative multiplication, which will be 
denoted by •. It induces a product on the tangent space of every point t e B. The 
class of the function F in O c defines a vector field E e T B , called the Euler vector 
field. 
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Given a holomorphic volume form to on (C 2/+1 , 0), possibly depending on t e B, 
we can equip p*O c with the so-called residue pairing: 



y 2m> J r d yo F---d y2l 



F 

where y = (yo, . . . , y2/) is a unimodular coordinate system for to (i.e. = dyo A • • • A 
Jy 2 /) and the integration cycle T e is supported on \d yo F\ = • • • = |5 V2; F| = e. Using 
that 7~g = p*O c , we get a non-degenerate complex bilinear form ( , ) on T B , which 
we still call residue pairing. 

For t e B and z e C, let be a semi-infinite cycle in C 2/+1 of the following 
type: 

S u 6 lim H 2M (£ 11+ \{RQ Z ~ l F(t,x) < -p};C) = C w . 

p— *oo 

The above homology groups form a vector bundle on B x C* equipped naturally 
with a Gauss-Manin connection, and S = ¥> t ,z may be viewed as a flat section. 
According to K. Saito's theory of primitive forms lfl"3l [T31 there exists a form o>, 
called primitive, such that the oscillatory integrals (d B is the de Rham differential 
on B) 



J s (t,z):=(2nzy l -hzd B ) f e^ FM toeT* B 



are horizontal sections for the following connection: 

v B,of = V^-z" 1 ^.,), l<i<N (10) 
Vfl /ai = d z -z- l e + z~ 2 E; . (11) 

Here V L c is the Levi-Civita connection associated with the residue pairing and 

e := V L - Q £ - (l - ^)ld, 

where d is some complex number. In particular, this means that the residue pairing 
and the multiplication • form a Frobenius structure on B of conformal dimension d 
with identity 1 and Euler vector field E. For the definition of a Frobenius structure 
we refer to [0 . 

Assume that a primitive form to is chosen. Note that the flatness of the Gauss- 
Manin connection implies that the residue pairing is flat. Denote by (r 1 , . . . , t n ) a 
coordinate system on B that is flat with respect to the residue metric, and write d t 
for the vector field d/dr'. We can further modify the flat coordinate system so that 
the Euler field is the sum of a constant and linear fields: 

N N 

E = Y J (l-ddT i d i + ]Tp l d i . 

i=l i=l 
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The constant part represents the class of / in H, and the spectrum of degrees 
d\ , . . . , d N ranges from to d. Note that in the fiat coordinates r ! the operator 
(called sometimes the Hodge grading operator) assumes diagonal form: 

e(d i ) = (^-d l )d i , l<i<N. 

Finally, the vectors v, e H appearing in formula © are the images of the flat vector 
fields dj via the Kodaira-Spencer isomorphism CQ). 

2.2. Period integrals. Given a middle homology class a e H 2 i{X G y, C), we denote 
by a U A its parallel transport to the Milnor fiber X t:A . Let d~ l a> be any 2/-form whose 
differential is co. We can integrate d~ x io over a t> \ and obtain multivalued functions of 
A and t ramified around the discriminant in S (over which the Milnor fibers become 
singular). 

Definition 2.2. To a e H 2 i(Xoy,C), we associate the period vectors lf(t,A) 6 
H (k e Z) defined by 

{If (t, A), d t ) := -{InT'dfdi f d-'oj , l<i<N. (12) 

Note that this definition is consistent with the operation of stabilization of sin- 
gularities. Namely, adding the squares of two new variables does not change the 
right-hand side, since it is offset by an extra differentiation (2^) -1 ^. In particular, 
this defines the period vector for a negative value of k > -I with / as large as one 
wishes. Note that, by definition, we have 

d^\t,A) = r* +l \t,X), keZ. 
The following lemma is due to A. Givental IfTOl . 
Lemma 2.3. The period vectors (fT2)) satisfy the differential equations 

<9,-/ w = -d i * t (d A I {k) ), l<i<N, (13) 

(A-E. t )d A I® = (e-k-^)l (k K (14) 

Using equation (fl4"l) . we analytically extend the period vectors to all \A\ > 6. It 
follows from (PT3l that the period vectors have the symmetry 

rf(t,A) = If(t-Al,0), (15) 

where t i-» t - Al denotes the time-A translation in the direction of the flat vector 
field 1 obtained from 1 6 H. (The latter represents identity elements for all the 
products • t .) 
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2.3. Stationary phase asymptotic. Let U0) (1 < i < N) be the critical values 
of F(t, •)• For a generic t, they form a local coordinate system on B in which the 
Frobenius multiplication and the residue pairing are diagonal. Namely, 

d/dui d/duj = djjd/diij , (d/duud/duj^ = Sy/Ai , 

where A, is the Hessian of F with respect to the volume form oj at the critical point 
corresponding to the critical value u t . Therefore, the Frobenius structure is semi- 
simple. 

We denote by % the following linear isomorphism 

: C N T t B , e t h> jArf/dut , 

where [e\, . . . , e N } is the standard basis for C^. Let U, be the diagonal matrix with 
entries u.\(i), . . .,u N (t). 

According to Givental O, the system of differential equations (cf. (fTOl) . CLD) 

zd i J(t,z) = d i * t J(t,z), \<i<N, (16) 

zd z J(t,z) = (8- z - 1 E. t )J(t,z) (17) 

has a unique formal asymptotic solution of the form ^fR^e 11 ' 11 , where 

R t (z) = l + R l (t)z + R 2 (t)z 2 + -- - , 

and Rk(t) are linear operators on C N uniquely determined from the differential equa- 
tions ([TBI and <fT71 ). Introduce the formal series 

f a (t,A;z) = Y J I a ) (t,A)(-z) k . (18) 

keZ 

Note that for A i -singularity F(t, x) = x 2 /2 + t we have u := u\(t) = t and the series 
(PT8T) takes the form 

f Al (t,A;z) = J]lf i (u,A)(-z) k , 

keZ 

where 

/*'(«, ,i) = ("I)' <2 y7|' 2 " ! U " «)" 4 "" 2 , t i 

^' V * = 2 (ST^! ( ' , -" )i *'' 2 - * i0 - 

The key lemma, which is due to Givental IfTOl is the following. 

Lemma 2.4. Let t e B be generic and B be a vanishing cycle vanishing over the 
point (t, Ui(t)) 6 Z. Then for all A near u t := u t {t), we have 

f f3 (t,A;z) = y t R t (z)e,f M (u h A;z). 
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3. Symplectic LOOP SPACE FORMALISM 

The goal of this section is to introduce Givental's quantization formalism (see 
[|9l) and use it to define the higher genus potentials in singularity theory. 

3.1. Symplectic structure and quantization. The space H := H((z -1 )) of formal 
Laurent series in z~ l with coefficients in H is equipped with the following symplec- 
tic form: 

Qtyufc) ■= Res^iC-z),^)) , 0i, 02 e H, 

where, as before, (, ) denotes the residue pairing on H and the formal residue Res z 
gives the coefficient in front of z~ l . 

Let {di}f =l and {<9 ! }^j be dual bases of H with respect to the residue pairing. Then 

Q.(d i (-zy k -\djZ l ) = 8 ij 8 kl . 

Hence, a Darboux coordinate system is provided by the linear functions q' k , ptj on 
*K given by: 

ql = n(d i (-zT k - 1 ,-), p Ki = Q(;d i z k ). 

In other words, 

oo TV oo N 

Hz) = J] Z + Z Z PkMSX-zT'- 1 , <pe<H. 

k=0 i=l k-0 i=\ 

The first of the above sums will be denoted (p(z)+ and the second <p(z)-. 
The quantization of linear functions on H is given by the rules: 

q[ = n- l/2 q[, Ki = ^ 2 ^j. 

dq\ 

Here and further, h is a formal variable. We will denote by Cr the field C((ft 1/2 )). 

Every <p(z) e fi gives rise to the linear function Q.((p, •) on "K, so we can define 
the quantization 0. Explicitly, 

oo N q oo N 

* = - %m Z Z ti&TT + r 1/2 Z Z puwd ■ (19) 

k=0 i=l k=Q ;=1 

The above formula makes sense also for 0(z) e H[[z, z -1 ]] if we interpret as a 
formal differential operator in the variables q' k with coefficients in Q,. 

Lemma 3.1. For all 0i , 02 e *H, we /zave [0i , 2 ] = ^X0i , 02)- 

Proof. It is enough to check this for the basis vectors d l (-z)~ k ~ l , diZ k , in which case 
it is true by definition. □ 
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It is known that the operator series %(z) '■= v F f i? f (z) v P~ 1 is a symplectic transfor- 
mation. Moreover, it has the form e A( - z) , where A(z) is an infinitesimal symplectic 
transformation. A linear operator A(z) on "K := //((z" 1 )) is infinitesimal symplectic 
if and only if the map (f> e *H i-» A(p e *H is a Hamiltonian vector field with a 
Hamiltonian given by the quadratic function h A ((p) = ^Q.(A(f), <p). By definition, the 

quantization of e A(z) is given by the differential operator e llA , where the quadratic 
Hamiltonians are quantized according to the following rules: 

_ d 2 • _ • d _ 1 • 

(PkjPij) =% —— (Pkjq J ,) =(q J ,Pk,i) =q J ,— , iff k q\) = -q\q\ . 
dq\dq\ oq k « 

3.2. The total ancestor potential. Let us make the following convention. Given a 
vector 

oo N 

q(z) = J] qkzk 6 H ^ ' qk = J] qi k di 6 H ' 

k=0 i=l 

its coefficients give rise to a vector sequence qo,qi, By definition, a formal 

function on H[z], defined in the formal neighborhood of a given point c(z) e H[z], 

is a formal power series in q - c ,qi - c\, Note that every operator acting 

on H[z] continuously in the appropriate formal sense induces an operator acting on 
formal functions. 

The Witten-Kontsevich tau-function is the following generating series: 

D pt (h;Q(z)) = exp(]T^- i f J(Q(^) + (20) 

g,n U - Jm s" (=1 

where Q , Qi, . . . are formal variables, and ^ (1 < i < n) are the first Chern classes 
of the cotangent line bundles on Ai gt „ (see |[T6l [T2l ). It is interpreted as a formal 
function of Q(z) = Zitto QkZ k 6 C[z], defined in the formal neighborhood of -z. In 
other words, £) pt is a formal power series in Q , <2i + 1 , Q 2 , Q3,. . . with coefficients 
in Cp)). 

Let t 6 5 be a semi-simple point, so that the critical values w,-(f) (1 < i < N) of 
F(t, •) form a coordinate system. Recall also the flat coordinates r = (r 1 ^), . . . , 7^(0) 
of t. The total ancestor potential of the singularity is defined as follows 

N 

q(z)) = ft, ] ] D pt (hAf, j q(z)) e q»[fe , <?i + 1, <? 2 . . . ]], (21) 
where ft,(z) := ^(z) 1 ^ 1 and 



00 N 



k=0 7=1 
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3.3. Proof of Theorem ll.il Using Lemma I2~4l it is easy to see that the ratio 

(/^(UXvM/^CUXl) 

is analytic in a neighborhood of A = uj for all m > 0. Furthermore, Y" J A %, = % Y u l , 
where Y u j- is the differential operator ([5]) on the variables Q J k := J 'qiJ -y/Kj (k > 0), 
defined for the A i -singularity F = x 2 /2 + Uj (see Lemma 6. 12 in Q). Using that the 
Virasoro operators in the Virasoro constraints for the point coincide with the polar 
part of YuJ (see Section 8.3 in [0), we get that Y" J A 3\ t {ti; q) is a formal series in q 
whose coefficients are analytic at A = Uj. This implies that the residue at A = Uj 
vanishes, which completes the proof. □ 

3.4. The Virasoro recursions. Let us compute the coefficient in front of q\ in 

The contribution from the j'-th residue is computed as follows. We put 
/3 = fij to avoid cumbersome notation. By definition the differential operator Y" J A is 
a sum of quadratic expressions of the following 3 types: 

(-If +k " <tf +1 \t, A), v a ) (/f +1) (;, A), v b ) Kd^, ' 

2(-l) r+1 (§ v \t, A), v a ) (lf' +l \t, A), v h ) q^d^ , (22) 

and 

(I ( /\t, A), v a ) (I ( /'\t, A), v b ) h~ l q a k , q\„ , 

where the sum is over all k',k" > and a,b = 1,2, ...,N. The only contribution 
could come from the terms (1221) . The coefficient in front of q\d q b is 

± (-l) fe+1 Res, =( , (t-- m -\t, A), vd (I« +1 \t, A), v b ). (23) 
Lemma 3.2. The following identity holds 

N 

£ Res, l=II , (lf\t, A), Va) (if'Xt, A), v b ) dA = 2(-l)*\»W,o, 

7=1 

for all k',k" e Z and a,b = 1,2, ... ,N. 
Proof. According to Lemma l2~4l we have 

oo 

lf(t,A) = Y j R l (-d z r l lf 1 (u j ,A)e j . 

1=0 

Using this identity we find that the j-th term in the above sum is 

oo 

CRi'Va, ej)( T R r v b , ej)(-l) l ' +l "Res A=Uj % +n (uj, A)lf^ l "\ Uj , A)dA. (24) 

/',/"=o 
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The above residue is non-zero only if k' + V = —k" - I". In the latter case using 
integration by parts we find that the residue is 

(-lf +l 'R e s^ Uj lf ) i (u j ,A)lf ) i (u j ,A)dA = 2(-lf +r . 

The sum (1241) becomes 



2(-lf J] CRi>v a ,ej)CRi»v b ,ej)(-lf6 y+l ^- k » 

l',t"=Q 

If we sum over all j = 1,2, . . .,N, since is an orthonormal basis of H, we get 

oo 

2(-lf Yj ( T Rl'Va, T Rl»V b )(-lf6 l , +l „,- k ,- k/ ,. 

rj"=o 

Using the symplectic condition R(z) T R(-z) = 1 we see that the only non-zero con- 
tribution in the above sum comes from the terms with /' = I" = 0, which completes 
the proof. □ 
The above Lemma implies that the coefficient <[2~3T ) is non-zero only if k = m and 
b = i and in the latter case the coefficient is 1 . In order to obtain a recursion relation 
for the correlators © we replace q\ = t\ - 1 and compare the genus g degree n 
(with respect to t) terms in the identity 

oo 

^(-l) m+1 (4" !+1) (r,^), v a )L m _ u ^;q) = 0. (25) 

m=0 

Note that if we ignore the dilaton shift, then (Y" ] A 3\,dA ■ dA)/Ji t (here • is the sym- 
metric product of differential forms) is a sum of terms of five different types. The 
first two are 

tjg-i 

— -(^'a,i ;i Ai),/ j a,i;(A2),t,...,t) , . (26) 

and 

Z , *M>*>...,t) , , +1 U+(t, A-MX-.., t) „ . ■ (27) 

t—i (n')'.(n")'. x I g',n'+\ \ / g",n"+l 

g'+g"=g 

n'+n"-n 

The other three types are 

Po ifi \t,A), (28) 

-—Q(/_ J (t,A;z),Kz))U^(t,A;^ l )X...,t) , (29) 

n ! N ' g>" 

and 

fr l a{(fj{t, A; z), t(z)) Q(/'(r, A; z), t(z)). (30) 
Let us point out that the ancestor correlators © are tame (see ifTOl ). which by 
definition means that they vanish if k\ + • • • + fc„ > 3g - 3 + n. In particular, the 
ancestor potential does not have non-zero correlators in the genus-0 unstable range 
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(g,n) = (0,0), (0, 1), and (0,2). However, motivated by the above formulas, it is 
convenient to extend the definition in the unstable range as well by setting 

(/'(U;<Ai),t) 02 := ft(/'(Uu),t(z)) 

(<f%(t, A; <f>{t, A- ^)} o>2 := P^ fi % A) 

and keeping the remaining unstable genus-0 correlators 0. If we allow such unstable 
correlators; then the terms (|29l ) and (1301) become the unstable part of the sum (1271) . 
while (1281) becomes the unstable correlator in the set (l26l) . 

The above discussion and the fact that the dilaton shift changes L m _! a simply by 
an additional differentiation -d/df;' n yields the following identities: 



+(*,//; 0i),t,...,t) 
1 N 



Res 



g,n+l 



A=Uj 



(I^\t,A),\)dA 



x 



7=1 



\ ' g-l,n+2 



g +8 =8 
n'+n"=n 



V r) (<f%(t, A; <Ai), t, . . . , t) (0j(f, A; <Ai), t, . . . , t) 

Z — i \W / e > n + 1 I g ,n"+l 



where (g, n) is assumed to be in the stable range, i.e., 2g - 2 + n > 0, we are 
allowing unstable correlators on the RHS, and we suppressed the dependence of 
the correlators on t e B ss . Note that the RHS involves differential forms that should 
be treated formally with respect to the symmetric product • of differential forms. 



dfi ■ dA ■ dA 
dA 



dfx ■ dA = dA - d/u.. 



The residue contracts dA, so at the end the RHS involves only djx. Let us point 
out that the tameness condition is crucial, because it implies that for the correlator 



insertion of the type ^ m Ig \t,A){-ij/) m , only finitely many terms contribute. In 



particular, although the infinite sum of differential operators in (1251) does not make 
sense in general, our argument goes through since on each step only finitely many 
terms of the sum (1231 ) contribute. 



(31) 
(32) 

(33) 



4. The Eynard-Orantin recursion 

The initial data for setting up the local Eynard-Orantin recursion is a complex 
line C with N marked points u\, . . . , u N and a certain set of 1- and 2-forms defined 
only locally. More precisely, for each i we have a multi-valued holomorphic 1-form 



EYNARD-ORANTIN RECURSION FOR THE ANCESTORS 15 

o)'(A) = P'(A) dA defined in some disk-neighborhood D, of u t , s.t., 

a-+i/2 



P\A) := Yj 



k=0 

For each pair we have a symmetric 2-form a)' j (ju, A) := P' j (/u,A)dA ■ dfi on 
D; x D 7 obeying the symmetry (z, //) <-> (7, A) and such that the function 

(// - ud m (A - Uj) i/2 P ij (M, A) = (m- ud m (A - Uj) i/2 P ji (A, M ) 

is analytic on D t x Dj except for a pole (with no residues) of order 2 along the 
diagonal in the case when i = j. In the latter case, we assume that the differentials 
are normalized in such a way that the Laurent series expansion with respect to fi in 
the annulus < \A - fu\ < \A - u\ has the form 

P U (M, A) = — + V Pl(A) ( M - Af. (34) 

Note that for each fixed k > 0, the functions P'l(A) are holomorphic on the punctured 
disk D* := Dj\ {w,} with a finite order pole at A = «j. 

Let us denote by X, the local system in a neighborhood of D, defined by the 
multi- valued function (A - Ui) l/2 . We define a system of symmetric multi- valued 
analytic differential forms af g ]n"' a "(A\, . . . ,A n ) for or* e £ ik and At e D* that are 
compatible with the (local) monodromy action on the local systems, i.e., the an- 
alytic continuation along a small loop around A^ = u ik transforms the differential 
form into aj^"^'"'""^, *n), where cr is the corresponding monodromy ac- 
tion on Xjj • 

Let a e £ t ,f3 e X, be any sections; then the base of the recursion is the following 

= 0, 

lpj(M,A)d/u-dA if(i,n)±(j,A), 
[Pl(A)dA-dA if = (j,A). 

The kernel of the recursion is the following ratio of 1 forms: 



(M, A) = 



I 4 P j (m, A ) dA du 

K^(m, A) = - ^ : ^, (35) 

^ ' 2 PJ(A) dA V 

where we fix /j. e D t \ {w,} and select a simple loop C,\ in Dj based at A that goes 
around Uj. Then the recursion takes the form 



co2f + l r a "(A , A u ...,A n ) = Y J Res,=„ ; K a ^(A , A) x 



./=! 
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?'+g"=g 

/'U/"={l,...,n} 

where we are assuming that 2g - 2 + n > 0, fij e _£y, the sum in the big brackets 
is over all splittings, n' and n" are the number of elements respectively in /' and 
/", and for a subset / c {1,2, . . .,n] we adopt the standard multi-index notation 
xi = (Xi t , . . . , Xi k ). Note that although the functions are multivalued, the local mon- 
odromy about A = uj leaves the expression invariant, so the residue is well defined. 

4.1. Proof of Theorem 11.21 In the settings of singularity theory, for a generic t e 
B ss we let the marked points be the critical values w, = Ui(t). The choice of a section 
of the local system X; is the same as choosing a vanishing cycle over A = w,. Let 
Wg)n '^ n (.Ai, . . . ,A n ) be the n-point series ©. 

In order to prove that these forms satisfy the Eynard-Orantin recursion, it is 
enough to notice that 

— {f + {t,A),tm =<fc(t,A;if,). 

Applying this formula n times to (T3TT) — (T33T) with (3 = /3i,...,J3„ we obtain the 
Eynard-Orantin recursion with 

co ij (ju, A) = [r + (t,»l ?-(t, A)], (36) 

and 

P j (A) = 4(I i p 1 \t,A),l), (37) 

where a, {3 are vanishing cycles vanishing respectively over /u = u t and A = Uj. Note 
that 

P ij '(M,A')dA' -tin = 2 [cT\(t,n),P-(t,A)] 



so the kernel is given by formula (1351) . 

In the opposite direction, in order to prove that the Eynard-Orantin recursion 
implies the Virasoro constraints, it is enough to notice that according to Lemma |3~2l 
we have the following identity 



1 N 

m = 2 X ReSi ="; Q(f ^' A > zX t(z)) A > ^ ■ D 



./=! 



Finally, let us point out that using Lemma l24l one can express the Laurent series 
expansions of co lj (p., A) and co J (A) in terms of the symplectic operator series %. The 
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answer is the following. Let W u e End(H) be defined via 

jrv MW y = '- r "- z > , 

k,l=0 ' 

then (// - w,) 1/2 (/l - «y) 1/2 P lj (pi, A) has the following Taylor's series expansion 



St 



Note that if i = j and we fix A near u{, then the Laurent series expansion of P'-'Qu, A) 
about fi = A does take the form (1341) . The Taylor's series expansion of is 

Up to an appropriate normalization of the correlation functions, our answer agrees 
with the formulas in 01. 
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